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Trapping cold atoms near carbon nanotubes: 
thermal spin flips and Casimir— Polder potential 

R. FermaniQ S. Scheel, and P.L. Knight 

Quantum Optics and Laser Science, Blackett Laboratory, Imperial College London, 
Prince Consort Road, London SW7 2AZ, United Kingdom 
(Dated: February 1, 2008) 

We investigate the possibility to trap ultracold 87 Rb atoms near the outside of a metallic carbon 
nanotube (CN) which we imagine to use as a miniaturized current-carrying wire. We calculate 
atomic spin flip lifetimes and compare the strength of the Casimir-Polder potential with the mag- 
netic trapping potential. Our analysis indicates that the Casimir-Polder force is the dominant loss 
mechanism and we compute the minimum distance to the carbon nanotube at which the atoms can 
be trapped to be larger than 100 nm. 
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I. INTRODUCTION 

Advances in magnetic trapping of clouds of ultracold 
atoms and Bose-Einstein condensates have received con- 
siderable research attention [ij, |2j, y, |4[ . The control and 
manipulation of atomic clouds is of fundamental impor- 
tance in the investigation of the basic physical properties 
of atom-surface interaction 0, d, 0K as well as in quan- 
tum information processing (QIP) Hi]. As trapped cold 
atoms appear to be very sensitive to magnetic-field vari- 
ations, they represent a powerful tool in magnetic field 
imaging as well as in gaining insight into atom-surface 
coupling phenomena [lO, EH- The challenge is to keep 
the atoms as close as possible to the substrate material 
to map the magnetic and electric fields in the vicinity of 
the surface. The combination of quantum state control 
with the development of ever-smaller magnetic traps is 
an essential element in the implementation of integrated 
quantum devices for fundamental research, quantum in- 
formation processing and precision measurement. 

Along with the push towards miniaturization evolved 
the idea of devising even smaller structures based on car- 
bon nanotubes (CNs) [HI]. Carbon nanotubes are carbon 
mono-layers rolled-up into cylinders of a few nanometers' 
diameter [l3|, [3]. They have been widely investigated 
theoretically and experimentally as they play a key role 
in miniaturized electronic, mechanical, electromechani- 
cal and scanning-probe devices. For their potential use 
as miniaturized current-carrying wires it is important to 
realize that the desired close proximity of a neutral atom 
and the carbon atoms that make up the nanotube can 
vastly enhance the influence of dispersion forces, and we 
address this point in the paper. 

It is well known that an atom held in a magnetic trap 
near an absorbing dielectric surface will be subject to 
thermally-induced spin flip transitions whose origin lies 
in fluctuating magnetic fields which can be attributed 



to resistive noise in the substrate [la, [lfj. In accor- 
dance with the fluctuation-dissipation theorem, dissipa- 
tion processes associated with a finite conductivity give 
rise to electromagnetic-field fluctuations. These fluctua- 
tions can be strong enough to drive spin transitions that 
lead to trapping losses [l?], El, Eil • 

In addition, an atom placed near a dielectric body will 
experience a dispersion force due to the presence of the 
dielectric material — the Casimir-Polder force [2(J Hi], 
[22l . [23j . The potential generating this force adds to the 
magnetic trapping potential and may cause the trap to 
become unstable at small distances. 

With this in mind, CNs seem to represent rather 
attractive structures for designing miniaturized mag- 
netic traps. This is on one hand due to the fact that 
they consist of a very small amount of dielectric mat- 
ter which means that unwanted dispersion forces such as 
the Casimir-Polder force are minimized. On the other 
hand, they also possess extremely homogeneous surfaces 
and are thus less likely to induce inhomogeneities in the 
potential surface of the trap. 

This paper is organized as follows. In Sec. [IT] we in- 
troduce the basic concepts of magnetic trapping of neu- 
tral atoms. In Sec. IHII the spin flip lifetime is calculated 
and compared to the tunneling lifetime resulting from 
the combination of the magnetic trapping potential and 
the Casimir-Polder potential. Both lifetimes are given 
in terms of the dyadic Green tensor. In order not to in- 
terrupt the flow of arguments and results we will present 
technical details about the conductivity and the Green 
tensor for a single wall CN in the Appendix. Some con- 
clusions are drawn in Sec. IIV1 



II. MAGNETIC TRAPS FOR NEUTRAL ATOMS 

In this section we briefly review magnetic trapping us- 
ing straightforward tools of electromagnetism. An atom 
with a magnetic dipole moment fj, placed in a magnetic 
field Bt experiences an interaction potential 
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V = -n-B T . 
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Assuming that a low-field seeking atom is in the state 
\IJFmp), Eq. ((T|) corresponds to a Zeeman energy V — 
gF/iB'nT'FB that depends on the quantum number mj? 
and on the magnitude of the field B = |Bx| while it is 
independent of the field's direction, with gp the Landc 
factor and fis the Bohr magneton. Let us consider a 
current / flowing through a wire along the z direction 
that generates a circular magnetic field B. 

Applying a homogeneous bias magnetic field Bf, point- 
ing in a direction orthogonal to the wire, a line of vanish- 
ing magnetic field parallel to the wire is created which is 
located a distance yo = (J-qI /(2nBb) away from the wire. 
The total field is given by Br = B + Bb whose gradient 
B'(r) — —fioI/(2irr 2 ), at the position of the field mini- 
mum ?/o, can be written as B'(yo) = b' = — 2B 2 ir/(noI). 
As B'(yo) is independent of position, the superposition of 
the magnetic field created by the wire and the homoge- 
neous bias field creates a two-dimensional quadrupolc- 
type trap (24J. In such a trap, the trapping poten- 
tial can be approximated by a linear function of the 
magnetic field gradient b' . The magnetic field can be 
expressed as B = b'xe x — b'ye y and its modulus is 
|B| = b'{x 2 + j/ 2 ) 1 / 2 = b'r with r denoting the distance 
from the trap centre. 

The modulus of the y-component of Bt vanishes at 
the trap center. In order to prevent Majorana transitions 
to non-trapped magnetic levels, a further offset field B„ 
parallel to the wire is applied with |B Q | <C |B|. The 
magnitude of the field near the center of the trap is then 
given by [25( 
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where the approximation holds for b'r <C B . The pres- 
ence of the B Q field changes the shape of the potential 
near the trap centre from being linear to harmonic. The 
interaction potential that follows from Eq. (TTJ is then, to 
a good approximation, a harmonic potential of the form 
V = Vq + \Muj 2 r 2 where M denotes the mass of the 
atom. The trap oscillation frequency Lu r is given by 



effects become important [27] . The physical properties 
of a CN are determined by the way in which the graphite 
sheet is rolled. The winding angle with respect to the 
hexagonal carbon lattice is usually described by two in- 
teger numbers (a, b) [H, H, HJ El . When 2a + b = 3n, 
where n is again an integer, a CN shows metallic be- 
haviour, otherwise it is semi-conducting. 

The axial conductivity <t zz (uj) and the resulting dielec- 
tric permittivity e(oj) of a (9, 0) carbon nanotube are 
calculated in Appendix [AJ At a frequency fo = 70 kHz, 
chosen as to correspond to an offset field B Q = 10~ 5 T = 
100 mG, we obtain <j zz (lu ) = 1.19 • 10 9 + 11.5i (Dm)" 1 
and e(ujo) ~ 3 • 10 14 i. Hence, a (9,0) carbon nanotube 
can indeed be considered as a metallic cylinder. 

Thermal fluctuations generate noise currents that lead 
to fluctuating fields near the body surface. We expect 
the noise due to these fluctuating fields to be much 
reduced in a CN compared to a dielectric bulk mate- 
rial due to the very small amount of matter involved. 
Nevertheless, thermal spin flips and the Casimir-Poldcr 
force cannot be neglected and needs investigation and 
a comparison of their effect. Both mechanisms origi- 
nate from the fluctuations of the electromagnetic field 
in the substrate. In particular, the spin-flip transitions 
are caused by the magnetic-field fluctuations while the 
Casimir-Polder force arises from both electric and mag- 
netic field fluctuations, the latter usually being negligible. 

In order to describe the two phenomena we utilize the 
quantization scheme of the electromagnetic field in the 
presence of dispersing and absorbing bodies [32j, [33j . As 
this theory is a macroscopic theory whose central quan- 
tities are linear susceptibilities, carbon nanotubes are 
probably at the limit to what we can actually describe 
with it. However, when viewed from distances that are 
several multiples of the bond lengths, the CN can be 
thought of as a homogeneous object so that the detailed 
structure from the surface cannot be resolved and QED 
in dielectrics can be safely used. This also assumes that 
the CN contains no impurities and shows no pitch alter- 
ations. 
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Moreover, B„ controls both splitting of the magnetic sub- 
levels by a frequency /o = \[iBB jh at the trap center 
as well as the stability of the resonance associated with 
the magnetic guide [2a |. 



III. TRAPPING LIFETIMES 

In this section, we investigate the two main limitations 
to the trapping lifetime, thermally-induced spin flip tran- 
sitions and the Casimir-Polder potential. In the follow- 
ing calculations, we take the current through the single- 
wall CN to be equal to / = 20 /iA which seems to be the 
largest current that can be sustained before saturation 



A. Spin flip lifetime 

If an atom is held sufficiently close to the CN surface it 
will experience quantum fluctuations of the electromag- 
netic field. At the center of the trap, the atom feels a 
constant magnetic field B Q . The atomic magnetic sub- 
levels are thus split by the Zeeman interaction and only 
a subset of these levels will experience an attractive force 
(low-field seeking states). An 87 Rb atom can be trapped 
in the hyperfine state \F, mp) = |2,2) is trapped, but 
only for sufficiently tight magnetic traps also in the state 
\F, tuf) — |2, 1). Transitions to lower magnetic sublevels 
allow the atom to escape. In the following, we disregard 
all the lower-lying states and treat the atom in the two- 
level approximation as the transition is the rate limiting 
step. 
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The lifetime of an atom due to spin flip transitions is 
given by the inverse of the spin flip rate [16j 
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where [Ib is the Bohr magneton, Sk is the kth vector 
component of the electronic spin operator, and g$ ~ 2 
the electron's g factor. Spin flips occur between the initial 
state \i) and the final state |/), the position r of the atom 
is taken to be the centre of the trap. 

The spin flip rate in Eq. ([4]) is given in terms of the 
dyadic Green tensor G(r, r, u>) which contains the physi- 
cal and geometrical information about the nanotube. We 
assume the CN to be in thermal equilibrium with the en- 
vironment at a temperature T. The total spin flip rate 
is then given by r to t = r(n t h + 1) where nth is the mean 
thermal occupation number nth = (e hu) °/ kBT — with 
ks denoting Boltzmann's constant. 
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FIG. 1: Spin flip lifetime of a rubidium atom near a (9,0) 
carbon nanotube with radius Rcn = 3.52 A. The trapping 
distance yo is varied between 1 and 200 nm. The other pa- 
rameters are: fo = 70 kHz and T — 380 K. 

In Fig. [T] we show the calculated spin flip lifetime 
tsf = 1/Ttot as a function of the trapping distance 
yo from the surface of a (9, 0) CN for a tempera- 
ture T = 380 K, corresponding to a thermal exci- 
tation energy of k B T = 5.2 ICT 21 J (~ 33 meV). 
We consider the ground state transition |2,2) — > |2, 1) 
for a 87 Rb atom with the transition frequency fo = 
ujq/2-k — 70 kHz. At such frequencies, the thermally- 
induced spin flips dominate the spontaneous spin flips 
as Hujo = 4.8 10~ 29 J (~ 0.3 neV) « k B T. We 
calculate the spin matrix elements relative to that tran- 
sition through the Clebsch-Gordon coefficients and ob- 
tain for the non- vanishing matrix elements |(i|Sa:|/)| = 

K*l4l/>l = i/4- 

To evaluate the Green tensor that satisfies the correct 
boundary conditions at the CN surface we used the for- 
mulas obtained in Appendix [Bl In particular, the electric 



surface current density creates a discontinuity in the tan- 
gential component of the magnetic field [34[ . To compute 
the full Green tensor, we use the method of scattering 
superposition of dyadic Green tensors that are expanded 
into cylindrical vector wave functions (see, e.g., [35l. |36I|) 
which somewhat differs from the approach employed in 

The lifetime increases with the atom-surface distance 
yo and follows the same power law encountered in [ijj] 
for a solid wire. According to Fig. [TJ at an atom-surface 
distance of approximately 20 nm a lifetime of the order 
of a few seconds is achievable. The spin flip lifetime can 
reach one minute for distances approaching 120 nm and 
exceeds more than 100 s for trapping distances larger 
than 160 nm. These results suggest that an atom can 
be held very close to a metallic CN for sufficiently long 
times, and this is in line with our expectations about spin 
flip occurrence and with the atom-loss rate estimations 
presented in [f2l |. 



B. Casimir— Polder potential 

The presence of macroscopic dielectric bodies changes 
drastically the structure of the vacuum electromagnetic 
field. One consequence is that an atom in its ground state 
placed sufficiently close to a dielectric body experiences 
a non-vanishing, in general attractive, dispersion force, 
the Casimir-Polder (CP) force [H El, IS E3|. Since 
the CP potential adds to the (repulsive) trapping po- 
tential, atoms can tunnel through the resulting potential 
barrier and get stuck at the nanotube surface. The life- 
time we have calculated in Sec. lIII Al provides information 
about the distance at which an atom can be held before 
thermally-driven spin flips occur in a given time, but the 
Casimir-Polder force may play an even bigger role for 
small enough distances. 

The Casimir-Polder potential can be derived in lowest- 
order perturbation theory within the framework of QED 
in dielectric media [39]. If we assume that an atom is in 
an energy eigenstate then the CP potential is given 
by the body-induced — i.e. dependent on the quantity 
of material — (and position-dependent) shift of the eigen- 
value AEi corresponding to this eigenstate \l). The CP 
potential can be expressed as [201 133 ] 



du u 2 af^ 1 



where iu = to and (to) is the atomic polarizability in 
lowest-order perturbation theory. In particular, for an 
atom in a spherically symmetric ground state, one finds 
that 

2 



(wz)Tr[G (s) (r,r, iu)} (5) 
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with djfc = (Z|d|fc) representing the matrix dipole el- 
ements relative to the transition from the atomic ini- 
tial state \l) to the allowed states with frequency 
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Uki = {Ef. — Ei)/h. The expression of the CP poten- 
tial in Eq. ([5]) is given in terms of the scattering part 
G( s \r, r, iu) of the Green tensor and the frequency in- 
tegral is performed along the imaginary axis. 

The Casimir-Polder potential has to be compared with 
the magnetic trapping potential in order to establish the 
size of its effect. In Fig. [2] we show Vtot, the total poten- 
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FIG. 2: The potential Vtot is the sum of the Casimir-Polder 
potential of Eq. <(Sj and the magnetic trapping potential of 
Eq. |T}. The plot represents Vtot for three different trapping 
distances yo = 100 nm, 150 nm, and 200 nm, respectively. 
A (9,0) CN is considered with a 20 /iA current, the spin-flip 
transition frequency is taken to be fo = 70 kHz. 

tial seen by the atom at three different trapping distances 
t/o = 100 nm, 150 nm, and 200 nm. Vtot is given by the 
the sum of the two potentials given in Eq. and Eq. ([5]). 
We assume a 20 fiA current is flowing through the CN 
as this seems to be the largest current a single-wall CN 
can potentially withstand [27] and the associated heating 
justifies our choice of the temperature to be significantly 
higher than room temperature. Among all the possible 
transitions \l) — •> \k), we only consider the lowest elec- 
tronic transition D2(5 2 Si/2 — » 5 2 P3/2) with wavelength 
A ~ 780 nm and dipole moment | d.2 1 = 4.227 ecto [ao: 
classical Bohr radius] . We assume that the D2 transition 
represents the main contribution to the atomic polariz- 
ability while others transitions bring about a negligible 
contribution to the CP force. 

In contrast to the spin-flip lifetime, tempera- 
ture effects are negligible here for two reasons. 
The resonant contributions corresponding to vir- 
tual dipole absorption [4l| are suppressed because 
hc/X = 2.5 10~ 19 J (~ 1.6 eV) > k B T. On the other 
hand, the spacing between the Matsubara frequencies is 
one order of magnitude smaller than the lowest electronic 
transition frequency (the relevant frequency range over 
which the CP potential has to be computed) so that the 
Matsubara sum can be replaced by the integral in Eq. |(5J). 

As it is evident from Fig. [2l Vtot forms a potential 
barrier whose height and width vary with the trapping 
distance yo. As mentioned previously, the addition of the 
offset field B Q changes the bottom of the potential well 



from linear to a harmonic trap which is, however, not 
visible on the scale of the figure. With decreasing yo the 
potential barrier becomes more and more shallow, until 
for atom-surface distances smaller than the critical value 
of yo ~ 100 nm the barrier effectively disappears. For 
trapping distances larger than that, the total potential 
shows a pronounced minimum. For example, for yo = 150 
nm we estimate the trap oscillation frequency to be u r ~ 
0.7 kHz, and the width and the height of the potential 
barrier to be 68.6 nm and 3.8 • 10~ 29 J, respectively. 
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FIG. 3: Tunneling time tcp as a function of the trapping 
distance yo, varied between 1 and 200 nm. 

Using the WKB approximation, we can estimate the 
tunneling probability T and the corresponding tunnel- 
ing lifetime tqp = 2-k / (Tuj r ). The result is shown in 
Fig. [3] for a ground-state atom trapped at varying dis- 
tances yo- From the comparison of Fig. [T] and Fig. [3l 
it is clear that the effect of the CP force cannot be ne- 
glected. For small enough atom-nanotube distances (and 
indeed for all distances shown in the figures) the tunnel- 
ing lifetime is several orders of magnitude smaller than 
the spin flip lifetime. For example, at a trapping distance 
yo = 150 nm we estimate tsf and tcp to be 94.4 s and 
0.2 s, respectively and a tunneling lifetime of a few sec- 
onds is achievable for trapping distances equal or bigger 
than 170 nm where the spin flip occurrence is no more a 
limiting factor. 



IV. CONCLUSION 

In this paper we have investigated a novel way of minia- 
turizing atomic magnetic traps by replacing solid current- 
carrying wires by carbon nanotubes as the elementary 
building blocks. At first sight, the advantages of using 
CNs are both their small diameter and the fact that they 
are effectively two-dimensional structures. Hence, one 
would expect from scaling arguments that traps at rather 
small atom-surface distances could be realized. 

We have investigated the loss mechanisms both due 
to thermally induced spin flips as well as by tunneling 
through the Casimir-Polder barrier. The calculations 
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have been performed within the framework of quantum 
electrodynamics in dielectric media which is valid as long 
as the dielectric properties of the nanotube can be de- 
scribed by a macroscopic permittivity, and if the experi- 
mental situation is such that the atomic structure of the 
CN cannot be resolved and macroscopic boundary con- 
ditions can be set. 

The spin flip lifetime has been found to scale accord- 
ing to our expectations. That is, this lifetime follows, as 
a function of the atom-surface distance, the same power 
law as in the case of a solid wire, with the result that 
for distances much larger than the radius of the nan- 
otube the expected lifetime exceeds several seconds. In 
contrast, the alterations of the trapping potential by the 
Casimir-Polder potential are much more severe. It ap- 
pears that the minimal feasible trapping distance is larger 
than 100 ran. The main reason explaining this result is 
that single-wall nanotubes cannot sustain high enough 
currents (and thus cannot generate deep enough mag- 
netic traps) as they saturate at high electric fields. As 
a potential remedy, it would be beneficial considering 
multi-wall nanotubes. An increased number of carbon 
layers would allow for higher current densities and con- 
sequently to a magnetic trapping potential that would 
be comparable with the Casimir-Polder potential even 
for smaller distances. 
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APPENDIX A: AXIAL CONDUCTIVITY AND 
DIELECTRIC PERMITTIVITY 

The calculation of the Casimir-Polder potential re- 
quires the full knowledge of the frequency dependence of 
the conductivity. In order to make our presentation self 
contained we quote some results on calculations of the 
axial conductivity that have previously been published 
elsewhere. 

Here we briefly review the frequency dependence of 
the axial surface conductivity a zz (uS) and of the dielec- 
tric permittivity e(u>) for a single- wall CN following the 
presentation in [29l . |30j|. For a dielectric medium, the 
linear relation between cr(u>) and s(u>) is s(u))/so — 1 = 
io{w) / {(jJEq), where e r (uj) — e(u>)/eo is the (complex) rel- 
ative dielectric permittivity. The Clausius-Mosotti equa- 
tion establishes the relation between the response of a 
medium to an applied field, i.e. the polarization, and its 
dielectric constant. Because of the cylindrical structure 



of CNs, their polarizability is highly anisotropic, with the 
principle axis of the polarizability tensor oriented parallel 
to the cylindrical axis 0, [H, H2, [H, HJ] . Consequently, 



one is allowed to neglect the azimuthal current [45| . The 
axial conductivity per unit length can be expressed as 

MM 



cr zz (R,u}) = - 



iujEQ e r (R, lo) — 1 
S p T 



(Al) 



where R = (Rcni </>, Z) is the radius vector of an arbi- 
trary point of the CN surface, S is the area of a single 
nanotube, pr is the tubule density in a bundle. 

The physical properties of a CN are determined by the 
way in which the graphite sheet is rolled. The winding 
angle with respect to the hexagonal carbon lattice is usu- 
ally described by two integer numbers (a, b). Depending 
on a and 6, CNs are either semi-conducting or metallic 
and particularly a CN exhibits metallic properties when 
2a + b = 3n, where n is again an integer [28|, H^, [30l. [3l|. 
An (a, b) CN has a one-dimensional bands 

E ± (N,p)=±t 



'2irN a + 26 , , , . 

x \ 1 1 + 4 cos ( — pi cos — + 4 cos^ — 

a 2a / 2 2 



. Vt , „ 2 Pi 



where £ is § times the interatomic distance, N = 
0, 1, . . . , a — 1, and iz/t < p < n/£ with p the wave num- 
ber. The corresponding Fermi distribution function is 
f(E) = l/{exp[/3(_E — (j,)] + 1} with inverse temperature 
j3 and chemical potential p. 

The main contribution to the conductivity is given 
by the dynamic conductivity due to the free carrier 
term e(.(oj) but for high frequency regimes another term 
£ b r {t-o), arising from the transition between the conduc- 
tion and the valence bands, becomes important such that 
the relative dielectric permittivity is given by e t (uj) — 
e((uj) + £ h r {uj) [2!j|45j]. The interband transition term is 
given by 



= 1 + 



eh 2 



at N 



*> l f(E + (N,p))-f(E-(N,p)) 
^ P E + (N,p) - E_(N,p) 

x [ReK (N,p)} 2 

(hw) 2 + ih 2 u)/T r - [E+(N,p) - E_(N,p)} 2{ ' 

where r r is a phenomenological relaxation time and pc — 
2ap T = (7rV|y(2i?cN^ 2 ) is the density of carbon atoms 
per volume [31] . The Drude term is given by 



huj(huj + iH/Tr) ' 

with LOpi the plasma frequency 



(A3) 



dp[lmK (N,p)f 



2_pc v r /e 

x{f(E+(N,p))+f'{E-(N,p))} 



(A4) 
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The quantity Kq(N,p) corresponds to the (dimension- 
less) matrix element of the momentum operator and 
is given in [2§|. The following parameters have been 
used in our calculations: to — 4.32 x 10~ 19 J, h/r r — 
4.8 x 1CT 21 J, I = 2.13 A, R CN = 3.52 A. 



where k = w/c and rj 2 = k 2 — h 2 . To enhance readability, 
we have omitted the tensor product symbol ® between 
the even and the odd cylindrical vector wave functions 
which are defined as 



APPENDIX B: GREEN TENSOR OF A 
SINGLE- WALL CARBON NANOTUBE 

In this section we present our calculation of the dyadic 
Green tensor for a single- wall CN. For a single- wall nan- 
otube, we can approximate the carbon layer by a bound- 
ary layer with zero thickness. In this way, the Green ten- 
sor exhibits a discontinuity in its first spatial derivative 
across the carbon layer. Due to its cylindric symmetry, 
the problem can be described adopting the cylindric ba- 
sis {e r , & v , e 2 } assuming the CN to be directed along e z . 
We use the method of scattering superposition (see, e.g. 
[35|,[36[). For an atom located in r' outside the CN, the 
Green tensor can thus be written as 

;) + G^ } (r,r» , r > R CN , 
,uj) ,r<R C N, 

(Bl) 

, to J is the unbounded (bulk) Green tensor 
representing the contribution of direct waves from the 
source at r' to the point r, and the two scattering con- 



G(r, 



where Go(r, r', 



G (r,r',u 
G^(r,r' 



tributions G^'(r, r', a;) and G^f ; (r, r', ui) describing the 
reflection and transmission of waves from/through the 
cylindrical surface. In order to satisfy the homogeneous 
Helmholtz equation and the radiation condition at infin- 
ity, the vacuum term and the two scattering terms can 
be taken to be in the following form [3(| 

G (r,r» = 

yt8(y — r') 
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The symbol Z n (x) has to be replaced either by the 
Bessel function J n (x) or, if the superscript (1) appears 
on the respective vector wave function, by the (outgoing) 

Hankel function of the first kind Hn\x). The primes 
in Eqs. (|B2j) - (|B4j) indicate the cylindrical coordinates 
(r' ,(f>' , z'). The coefficients C m p (to = 1, 2, 3 and 4, 
and P = H , V) need to be determined from the bound- 
ary conditions for the electric and magnetic field compo- 
nents on the CN surface. The electric field satisfies the 
boundary condition 



e r x E(r,cj)| 



E(r,w) 



r=Rr 



= 0, (B7) 



while the electric surface current density creates a dis- 
continuity in the tangential component of the magnetic 
field 



e r x H(r,w)| + -H(r,u;)| = J(r, u)\ r = Rat} . 

(B8) 

Equations (|B7[) - (|B8[) translate into the respective bound- 
ary conditions for the Green tensor 



G(r,r»| r=fl± __-G(r,r»| r=B __ =0, (B9) 



r=R- 



r xVx G(r,r',oj)\ D+ -G(r,r',w)| „_ 

iu>fi cT(r)-G(r,r',u)\ r=Rc $10) 



where er(r) is the (diagonal) conductivity tensor whose 
only nonzero element is a zz (R, to). 



Substituting the decomposition (|B 1[) . together with 
Eqs. (|B2|) - (IB4[) . into the boundary conditions (|B9[) and 
(jBlpp leads to two sets of four equations for each po- 
larization H and V, that enable us to determine the 16 
coefficients C m p, 



7 



2 2 

- —H n (nr)C 2H + -rJ n {w)CiH = 0, 

thfi ihfi 
-d r H n (rjr)CiH ± -^H n (rjr)C 2 H +9 r J n (rir)C3HT — J n {w)C±H = d r J n (nr), 
kr kr 



-n H n (nr)C 1H + V J n {w)CsH 
T H n (r)r)CiH - kd r H n {nr)C 2H ± J n {w) C 3H 



+ ( kd r J n {nr) - iu)^ Q a zz —J n {r}r) I C. 



-4H 



r] 2 J n (r)r), 



ihn 

= ± JnW") 

r 



(Bll) 

(B12) 
(B13) 

(B14) 



and 



—H n (nr)C lv 
k 



ihn 
— H n (r]r)Civ 
kr 



ihn ihn 
d r H n (nr)C 2V ± -, — Jniw^v + d r J n {nr)C iV = ±-j — Jniw), 
kr kr 



ihn I rj~ 

-kd r H n (nr)C lv ± —H n (nr)C 2V + ( kd r J n (r]r) - iujn a zz — J n (nr) ) C 3V 



-n 2 H n (nr)C 2V + n 2 J n (nr)C iV 

t 
k 

ihn 

T J n \J]r)C w 

r 



0, 



kd r J n (nr). 



(B15) 

(B16) 
(B17) 



(B18) 



r 



The appearance of the axial conductivity a zz (R, to) in the 
boundary conditions (|B14[) and (|B18|) reflect the jump 
condition (|B10|) of the derivative of the Green tensor at 
the boundary layer. 

As we need to compute only the Green function in 
the region where the atom is located, only C( 12 )p need 
to be determined. On using various properties of the 
Bessel functions such as the Wronskian between the 
Bessel function J n (x) and the Hankel function H^\x), 

J n (x)Hi iy (x) - j;(x)iji 1} (a;) = 2/(irx), we obtain that 
the only non-zero coefficient is 



Finally, the Green tensor for an atom located at a posi- 
tion r' outside the CN can be expressed as 



G(r,r',w)=G (r,r» 



n=0 



(B20) 



Civ — — 



Trn LoRcNcr zz ri 2 J 2 (nRcN) 



2k 2 + nn u;R CN a zz ri 2 J n (nRcN)H n (nRcN) 



(B19) 



Equation (|B20[) . together with Eq. (|B2p . is the expression 
for the Green tensor used throughout this article. 
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